Dark torsion as the cosmic speed-up 
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It is shown that the recently detected acceleration of the universe can be understood by considering 
a modification of the teleparallel equivalent of General Relativity (TEGR), with no need of dark 
energy. The solution also exhibits phases dominated by matter and radiation as expected in the 
standard cosmological evolution. We perform a joint analysis with measurements of the most recent 
type la supernovae (SNe la), Baryon Acoustic Oscillation (BAO) peak and estimates of the CMB 
shift parameter data to constraint the only new parameter this theory has. 



I. INTRODUCTION 



The discovery of an unexpected diminution in the ob- 
served energy fluxes coming from type la supernovae 
P, [2] has opened one of the most puzzling and deep- 
est problems in cosmology today. These observations 
have been interpreted as solid evidence for an acceler- 
ating universe dominated by something called dark en- 
ergy. Although the cosmological constant seems to be 
the simplest explanation for the phenomenon, several 
dynamical scenarios have been tried out since 1998 (see 
e-g-, [1, 0) Hi)- While some authors sustain the idea of 
the existence of a dark energy, others propose modifica- 
tions of the Einstein-Hilbert Lagrangian known as f{R) 
(11, 0, S i 0, [H [H [H, H or [il, M [13 for recent 
reviews) as a way to obtain a late accelerating expan- 
sion. A great difficulty these theories have, from the 
point of view of the metric formalism, is that the re- 
sulting field equations are 4th order equations, which in 
many cases makes these hard to analyze. Besides, the 
simplest cases of the kind f{R) = R — BJR^ have shown 
difhculties with, weak field tests [H, gravitational 
instabilities [20| and do not present a matter dominated 
era previous to the acceleration era [U, [S^. Alterna- 
tively, the Palatini variational approach for such f{R) 
theories leads to 2nd order field equations, and some au- 
thors have achieved to put observational constraints to 
these theories [H, [23, H^- However in many cases the 
equations are still hard to work with, as evidenced by 
the functional form of the modified Friedmann equation 
for a generic f{R)- Recently, models based on modified 
teleparallel gravity were presented as an alternative to 
inflationary models [26l . [27j . In this paper we show a cos- 
mological solution for the acceleration of the universe by 
means of a sort of theories of modified gravity, namely 
/{Lt), based on a modification of the teleparallel equiv- 
alent of General Relativity (TEGR) Lagrangian j28l . [29j 
where the torsion will be the responsible of the observed 
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acceleration of the universe, and the field equations will 
always be 2nd order equations. 



II. GENERAL CONSIDERATIONS. FIELD 
EQUATIONS 

Teleparallelism [2^,[2l] uses as dynamical object a vier- 
bein field ei{x^), i = 0,1,2,3, which is an orthonormal 
basis for the tangent space at each point of the man- 
ifold: ■ Bj — rjij, where rjij = diag{l, —1, —1, 
Each vector can be described by its components ef , 



fi = 0,1,2,3 in a coordinate basis; i.e. 



Notice that latin indexes refer to the tangent space, 
while greek indexes label coordinates on the manifold. 
The metric tensor is obtained from the dual vierbein as 
9fj.:^ix) = rji j e''^{x) el{x) . Differing from General Rela- 
tivity, which uses the torsionless Levi-Civita connection, 
Teleparallelism uses the curvatureless Weitzenbock con- 
nection [301, whose non-mdl torsion is 



T =T - r = e- 

/il^ Z^/i fllf ^? 



(1) 



This tensor encompasses all the information about the 
gravitational field. The TEGR Lagrangian is built with 
the torsion ([T]) , and its dynamical equations for the vier- 
bein imply the Einstein equations for the metric. The 
teleparallel Lagrangian is [H, [3l|, [13] , 



where: 



and iiT^p is the contorsion tensor: 
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(2) 



(3) 



(4) 



which equals the difference between Weitzenbock and 
Levi-Civita connections. 

In this work the gravitational field will be driven by a 
Lagrangian density that is a function of Lt- Thus the 
action reads 



1 
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d^x e /(Lt) 



(5) 
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where e = det{e^^) = ^J~^. The case /{Lt) = Lt cor- 
responds to TEGR. If matter couples to the metric in 
the standard form then the variation of the action with 
respect to the vierbein leads to the equations 



1 
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+ 5/''9,(LT)/"(iT) + 7</(^T) = 



(6) 



where a prime denotes differentiation with respect to Lt, 
S^'^^ — e/Sp'^'^ and T^^ is the matter energy- momentum 
tensor. The fact that equations ([6]) are 2nd order makes 
them simpler than the dynamical equations resulting in 
f{R) theories. 



III. COSMOLOGICAL SOLUTION AND 
OBSERVATIONAL CONSTRAINTS 

We will assume a flat homogeneous and isotropic FRW 
universe, so 



e* = diag{l, a{t), a{t), a{t)) 



(7) 



where a{t) is the cosmological scale factor. By replacing 
in (d]), (01 and (g]) one obtains 



Lt = SP^^Tp., = -6 



(8) 



H being the Hubble parameter H ^ da^^. As a remark- 
able feature, the scale factor enters the invariant Lt just 
through the Hubble parameter. The substitution of the 
vierbein ([71) in © for i = = yields 

12H^ f'{LT) + f{LT)^lQTTGp (9) 

Besides, the equation i = 1 = is 

A8H^f"{LT)H - f'{LT)[l2H^ + AH] - f{LT) = 167rGp 

(10) 

In Eqs. (pIfTOl). p{t) and p{t) are the total density and 
pressure respectively. It can be easily derived that they 
accomplish the conservation equation 



dt 



(11) 



whatever /{Lt) is. Thus, if the state equation is p = wp 
then p evolves as p a (1 -f (z is the cosmological 

redshift). 

We are interested in obtaining an accelerated expan- 
sion without dark energy but driven by torsion. For this 
we will try with a kind of /{Lt) theories: 



/{Lt) = Lt 



{-Lrr 



(12) 



being a and n real constants to be determined by ob- 
servational constraints. Although the functional form of 



p2)) is similar to those considered in /{R) literature, now 
the guideline towards modified gravity is H instead of R. 
This fact gives to these theories another interesting fea- 
ture because H is the most important cosmological vari- 
able. For later times the term —a/{-~LT)'^ is dominant, 
while in early times, when — s- oo, General Relativ- 
ity is recovered. From ^ along with the modified 
Friedmann equation results to be 



(2n-t- 1) a _ 8 



irGp 



(13) 



(a functional dependence similar to the results other au- 
thors arrived, through different theoretical motivations 
such as m, Q ) . 

Now, replacing p — Pmo(l + z)^ + Pro(l + z)'^, and 
calling Qi = SttG pio/ {3H^) the contributions of matter 
and radiation to the total energy density today, Eq. ([13]) 
becomes 



y-{y-B) = C 



(14) 



where y = H^/H^, B = f7„,(l + z)^ + f7,.(l + z)^ and 
C = a(2n+ 1) (6i/2)-("+i). The evaluation of this equa- 
tion for 2; = allows to rephrase the constant C as a 
function of f2i and n: C = 1 — ~ ^r- For a = (then 
1 = Vim + ^r) the GR spatially flat Friedmann equation 
= H^B is retrieved. The case n = recovers the GR 
dynamics with cosmological constant £7 a = 1 — — ^r- 
Notice the functional simplicity of (fT3| compared with its 
analogue in /{R) theories. Compared with GR, n is the 
sole new free parameter in (jl4p , since specifying the value 
of n and {^r) the value of a (in units of jg 
automatically fixed through the relation (fT5|) . In order 
to obtain H{z) we solve numerically the equation (fT4|) . 

Since the most solid evidence for the acceleration of 
the universe comes from measurements of luminosity dis- 
tances for type la supernovae, we will use the most recent 
compilation of 307 SNe la events (the Union sample) [2| , 
to put constraints in the n— plane. The predicted dis- 
tance modulus for a supernova at redshift for a given 
set of parameters P=(n, $!,„), is 



p(z I P) 



M = blog{dL) + 2b 



(15) 



where m and M are the apparent and absolute magni- 
tudes respectively, and d^ stands for the luminosity dis- 
tance (in units of megaparsecs) , 



dL(z;P) = {l + z) f 
Jo 



dz' 



H{z',P) 



(16) 



where H{z; P) is given by the numerical solution of (jl4p . 
We use a statistic to find the best-fit for a set of pa- 
rameters P (marginalizing over Hg), 



XsNe 



Mz\P)-^,f^{z) 



(17) 



where pi{z \ P) is defined by (fT^ . p°''^ and cr^ are the 
distance modulus and its uncertainty for each observed 
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value i. As it is known, the measurements of SNe la 
are not enough to constraint thoroughly. To perform 
the statistic we also consider, on one hand, the informa- 
tion coming from the BAO (Baryon Acoustic Oscillation) 
peak detected in the correlation function of luminous 
red galaxies (LRG) in the Sloan Digital Sky Survey [ssl ]. 
The observed scale of the peak effectively constraints the 
quantity (assumed a AC DM model), 



^0.35 = i?v(0.35) = 0.469 ± 0.017 (18) 



where z = 0.35 is the typical redshift of the LRG and 
Dv is defined as 



Dv{z) 



dz' \2 



-,1/3 



H{z) \Jo H{z') 



(19) 



On the other hand, we have also included in the statis- 
tic the CMB shift parameter, which relates the angular 
diameter distance to the last scattering surface with the 
angular scale of the first acoustic peak in the CMB power 
spectrum. In order to do this, we have considered a radia- 
tion component Vlr—b x 10^^. The CMB shift parameter 
is given by [36| . 



/.1089 

^1089 — y^hnHl / 
Jo 



dz 
H{z) 



= 1.710 ±0.019 (20) 



We can use both parameters since our model presents 
matter domination at the decoupling time. Figure [1] 
shows the Hubble diagram for the 307 SNe la belong- 
ing to the Union sample. The curves represent models 
with values of Q,m and n obtained from minimizing us- 
ing only SNIa and SNIa+BAO+CMB as well. It is also 
shown as a reference the ACDM model with fi^ — 0-26. 
The obtained values for the best-fit to the SNe la data 
only are fim = 0.42 and n = 1.30 with the reduced 
xl = xLiJ" - 1-02 (or equivalently, Ax^„ = -1.1) 
where i> is the number of degrees of freedom. 

In Figure[2]wc plot the distance modulus residual (A/i) 
from the ACDM model to better appreciate the discrep- 
ancies between our model and ACDM. 

Figure [3] shows the confidence intervals at 68.3%, 
95.4% and 99.7% for the joint probability of the param- 
eters n and f2„i, having combined the SNe la data with 
BAO and CMB parameters. This analysis yields that 
the best-fit to all data is achieved with n = —0.10 and 
a„ = 0.27 (with a xLn/'^ - 1-01, AxL„ = -1.2) and 
also the values of the parameters lie in the ranges (at 
68.4% C.I.): n £ [-0.23, 0.03]; a« e[0.25, 0.29]. 

For our model we have analyzed as well the total and 
effective equations of state as a function of z. From IT^l) 
and ([TU)) along with , one can define a torsion contri- 
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FIG. 1: Hubble diagram for 307 SNe la from the Union sam- 
ple 01 • The curves correspond to the concordance model 
ACDM with nA=0.74 and f2m=0.26 (dotted line), and our 
models with the values corresponding to the best-fit r2m=0.42 
and n=1.30 (dashed line) and also with those coming from the 
joint analysis of SNIa+BAO+CMB, Q^=0.27 and n=-0.10 
(solid line). 



^ -0.10. 



best-fit SNIa 

best-fit SNIa+BAO+CMB 



FIG. 2: Distance modulus residual from the ACDM model 
for the same values from Figure [1] 



bution to the density and pressure as 
3 (2n+l)a 



Pt 
Pt 



a 



{6H^)-('''+^) H[4n{n + 1) - 2n] 
(6iJ2)-"- 



87rG 

- 6n(6i/2)-(«+i)i/2 _ 

to rewrite the dynamical equations as 
iirC, 



i/2 = 



a 
a 



-(p + Pt) 



Then, by using 



SttG 
and 



[p + PT + 3{p + pt) 



(21) 

(22) 
(23) 



the total and effective equa- 
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FIG. 3: Confidence intervals at 68.3%, 95.4% y 99.7% in the 
n — Qrn plane coming from combining SNe la, BAO and CMB 
data. The best-fit to this joint analysis is reached with the 
values n — —0.10 and Sim = 0.27. 



tions of state are written as, 
P + Pt 



Wtot 



Pt 



2{l + z) dH 
'3H lb 



w. 



eff 



PT 

Pt 



(24) 



(25) 



Figure S] shows the evolution of the total equation of 
state Wtot as a function of z for our model with the val- 
ues of the best-fit. There can be observed the last three 
phases of the evolution of the universe: radiation dom- 
inated {w — 1/3), matter dominated {w = 0) and late 
acceleration {w ~ —1). Figure [5] shows the effective 
equation of state coming from the dark torsion contribu- 
tion. Finally, analyzing we find that our model with 
n = —0.10 predicts that the transition from deceleration 
to acceleration occurs at Zacc — 0.74 in good agreement 
with recent works [s^]. 

An interesting point to be highlighted is that equation 
(fT5)l reveals that a value of n > 0, as the one obtained 
by considering only SNIa data, implies that the effec- 
tive dark torsion is of the phantom type [11] . That is, 
since H decreases toward the present time, the dark tor- 
sion density increases instead of diluting with expansion 
(weff < —1). However, when combining the complete 
data with SNIa+BAO+CMB we can see from Fig. [3] 
that it is slightly favored (Icr c.l.) a model with n ^ 0. 



IV. CONCLUSIONS 

A theory /(ir) based on a modification of the telepar- 
allel equivalent of General Relativity (TEGR) -where tor- 
sion is the geometric object describing gravity instead of 
curvature and its equations are always of 2nd order- is 
remarkably simpler than f{R) theories. We have tested 



FIG. 4: The curves correspond to the total equation of state as 
a function of z expected for the standard concordance model 
ACDM with f2A=0.74 and f2m=0.26 (dashed line), and for 
our model (solid line) with the values of the best-fit coming 
from SNIa-HBAO-hCMB, ^^=0.27 and n=-0.10. Three cos- 
mological phases observed. 
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FIG. 5: Effective equation of state as a function of z 
for our model with the values of the best-fit coming from 
SNIa+BAO+CMB, Q^=0.27 and n=-0.10. 



the theory /{Lt) = Lt — q;(— Lt)^" with the aim of re- 
producing the recently detected acceleration of the uni- 
verse without resorting to dark energy. We have here 
performed observational viability tests for this theory by 
using the most recent SN la data, and combined them 
with the information coming from BAO peak and CMB 
shift parameter in order to find constraints in the n — flm 
plane. At 68.3% c.l. we found that the values lie in the 
ranges n G [-0.23,0.03] and e [0.25,0.29]. The val- 
ues for rim are consistent with recent estimations ob- 
tained by other authors (see, e.g., ^). The model 
with the best-fit values minimizing the that combines 
SNIa+BAO+CMB data (n = -0.10 and fl^ ^ 0.27) 
exhibits the last three phases of cosmological evolution: 
radiation era, matter era and late acceleration; this last 
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stage having started at Za 



0.74. 
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